reg_data_10mi_raw <-
here("Shared/Data/WaterAnalysis/reg_data_10mi_raw.rds") %>%
readRDS()
# - Data on TB - #
TB_10mi <- reg_data_10mi_raw[nrdname=="Tri-Basin",]
# - Data on LR - #
LR_10mi <- reg_data_10mi_raw[nrdname=="Lower Republican",]
nrd_owner_name and
year.acres.LR_10mi_By_y_owner <-
LR_10mi %>%
.[,.(
usage = weighted.mean(usage, acres/sum(acres)),
treat2 = mean(treat2),
# --- soil --- #
silttotal_r = weighted.mean(silttotal_r, acres/sum(acres)),
claytotal_r = weighted.mean(claytotal_r, acres/sum(acres)),
slope_r = weighted.mean(slope_r, acres/sum(acres)),
ksat_r = weighted.mean(ksat_r, acres/sum(acres)),
awc_r = weighted.mean(awc_r, acres/sum(acres)),
# --- weather --- #
pr_in = weighted.mean(pr_in, acres/sum(acres)),
pet_in = weighted.mean(pet_in, acres/sum(acres)),
gdd_in = weighted.mean(gdd_in, acres/sum(acres)),
# --- tr --- #
cntyname_fix = unique(cntyname_fix)
), by = .(nrd_owner_name, year)] %>%
.[,county_year := paste0(cntyname_fix, "_", year)]
se_vars <- names(LR_10mi_By_y_owner)
TB_10mi_By_y_well <-
TB_10mi %>%
.[,county_year := paste0(cntyname, "_", year)] %>%
.[, ..se_vars]
reg_data_10mi <- bind_rows(LR_10mi_By_y_owner, TB_10mi_By_y_well)
tr=="5_22"
tr=="5_22"# === LR data === #
LR_10mi_By_owner <-
LR_10mi_By_y_owner %>%
.[,.(
usage = mean(usage),
treat2 = mean(treat2),
# --- soil --- #
silttotal_r = mean(silttotal_r),
claytotal_r = mean(claytotal_r),
slope_r = mean(slope_r),
ksat_r = mean(ksat_r),
awc_r = mean(awc_r),
# --- weather --- #
pr_in = mean(pr_in),
pet_in = mean(pet_in),
gdd_in = mean(gdd_in),
# --- tr --- #
cntyname = unique(cntyname_fix)
), by = .(nrd_owner_name)] %>%
.[,nrd_owner_name := NULL]
# === TB data === #
TB_10mi_By_well <-
TB_10mi %>%
# remove data points related to tr=="5_22"
.[tr!="5_22",] %>%
.[,.(
usage = mean(sum(volaf*12)/sum(acres)),
treat2 = mean(treat2),
# --- soil --- #
silttotal_r = mean(silttotal_r),
claytotal_r = mean(claytotal_r),
slope_r = mean(slope_r),
ksat_r = mean(ksat_r),
awc_r = mean(awc_r),
# --- weather --- #
pr_in = mean(pr_in),
pet_in = mean(pet_in),
gdd_in = mean(gdd_in),
# --- tr --- #
cntyname = unique(cntyname)
), by = .(wellid)] %>%
.[,wellid := NULL]
# unique(TB_10mi_By_well$treat2)
agg_reg_data_10mi <- rbind(LR_10mi_By_owner, TB_10mi_By_well)